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Abstract 



We include two loop, relativistic one loop and second order relativistic tree level 
corrections, plus leading nonperturbative contributions, to obtain a calculation of 
the lower states in the heavy quarkonium spectrum correct up to, and including, 
0(a 4 ) and leading A 4 /m 4 terms. This allows us, in particular, to obtain a model 
independent determination of the pole masses of the b, c quarks, 

m b = 5 0151 to MeV; m c = 1 884±IH MeV 

to which correspond the MS masses, 

fh b (ml) = 4453±|2 MeV; m c (ml) = 1 547t^jj MeV. 

The decay _T(T — > e + e~) is found in agreement with experiment, 

r(T -> e+O = 1. 135^-jg keV (exp. = 1.320 ± 0.04 keV), 

and the hyperfine splitting is predicted to be 

M(T) - M(rf) = 48.5tl2:2 MeV . 
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1 Introduction 

In recent years it has become possible to perform rigorous QCD analyses of heavy quarkonium systems, and 
this due to two reasons. First, radiative corrections have been calculated to increasing order of accuracy. 
The one loop corrections to the nonrelativistic (NR) spin-independent (SI) potential were calculated already 
in 1980 [1] . This was extended in refs. 2, 3 to the spin-dependent corrections, and in ref. 4 by including the 
velocity corrections to the SI part. Finally, the two loop nonrelativistic, spin independent corrections to the 
potential have been evaluated rece ntly[ 5 L 

Secondly, Leutwyler and Voloshin (ref. 6; see also ref. 7) have shown how to take into account, to leading 
order, nonperturbative (NP) effects, associated with the nonzero value of various condensates, of which the 
leading contribution is that of the gluon condensate (a s G 2 ). This has been implemented, together with the 
potential obtained with radiative corrections to one loop, in refs. 4, 8 where a study of bound states bb with 
nl = 10, 20, 21 and cc sates with nl = 10 was given 1 . The analysis was extended in refs. 9, 7 with the 
inclusion of size effects and higher condensates. 

The overall conclusion of these analyses is that pure QCD, without recourse to introducing phenomeno- 
logical interactions, produces a good description with manageable errors of the bb ground state and, to a 
lesser extent, of the splitting M(T) — M{j]b) and the decay T — > e + e~ . The description of the ground state of 
cc and of the excited states n = 2, I = 0, 1 of bb was shown to be even less reliable: the corrections are large, 
in some cases larger than the nominally leading terms. Still, it was possible, by using the renormalization 
point (i as a free parameter, to get a fairly accurate description of all n = 2 states including tensor and LS 
splittings! 8 ] _ 

In the present paper we extend this analysis by including the two loop corrections to the SI, NR potential 
recently calculated^ adding also velocity corrections to certain one loop pieces to get a calculation accurate 
up to, and including, 0(a 4 ) corrections. By taking into account the leading nonperturbative terms, we also 
include in the analysis corrections of order A 4 /m 4 . The main improvement so obtained is that we get an 
extremely stable and precise determination of the ground state of bb and, to a lesser extent, cc. If we invert 
the calculations we can deduce quark masses from the masses of the X, J ftp particles. We then find, for the 
pole masses, 

m t = 5 015l^° MeV, 

m c = 1884^22 MeV ^ 

m b (m 2 b )= 4453^° MeV, 
fh c (m 2 c ) = 154711$ MeV. 

The error includes the estimated theoretical error of the calculation, see the text. Note that (1.1) are very 
precise: as stated they are correct to order a 4 , and leading, 0(A 4 /m 4 ) nonperturbative effects. This is to 
be compared with estimates based on sum rules! 10 ! which are only accurate to order a 2 , or previous bound 
state calculations^, accurate only to third order in a s . 

The decay r(T — ► e + e~) is also given in agreement with experiment, within errors: 

r(Y -> e+e") = 1.135lg;^ kcV (exp. = 1.320 ± 0.04 kcV). 

The hyperfine splitting is predicted to be 

M(T) - M(n) = 48.5±J|;2 MeV . 

For higher states (nl = 20, 21) the errors are much larger but, within these, one has compatibility with 
experiment (cf. Subsect. 4.3). 

2 The Effective Potential 

We follow the method of effective potentials of Gupta et al' 3 ', and the renormalization scheme of ref. 4. The 
Hamiltonian for quarkonium may then be split in the following form: 

H = H<®+Hi (2.1a) 



to which correspond the MS masses, 



1 Throughout this paper n, I will denote the principal quantum number and the angular momentum of bound states. 
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where is 



a s (n 2 ) =a s (fi 2 ) \ 1 + ( ai + 



+ 



2 / 7T 
^2 



and can (and will) be solved exactly. i?i may be written as 

Hi = Free + VI L) + F 2 (i) + + V s . rel + F pin , 



(2.1b) 



(2.1c) 
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(2.1d) 



Here the running coupling constant has to be taken to three loops. For the values of the constants entering 
above formulas, cf. the Appendix. a\ was calculated in ref. 1, ai in ref. 4 and b\ in ref. 5. The other terms 
in (2. Id) can be obtained by use of the renormalization group, see e.g. ref. 4, or are well-known tree level 
relativistic corrections (including kinetic energy corrections). 

A few words are due on Eqs. (2.1). First of all, they only take into account the perturbative part of the 
calculation; NP effects will be incorporated later. Secondly, it should be noted that Hi contains a velocity- 
dependent one loop piece, F.rci- This is because the average velocity in a Coulombic potential is (|t>|) ~ a s , 
hence a calculation correct to order requires tree level 0(v 2 ) and one loop 0(H) contributions. All 
these terms in Hi may be treated as perturbations to first order, except V^ L \ For this, the second order 
perturbative contribution is required as this also produces a correction of order a*. 

A last comment concerns the renormalization scheme. We have followed ref. 4 in renormalizing a s in the 
MS scheme; but the mass m that appears in Eqs. (2.1) is the two loop pole mass. That is to say, it is defined 
by the equation, 

S^{j) = m,m) = (2.4) 

where S2 {i> ,rn) is the quark propagator to two loops. One can relate m to the MS parameter, also to two 
loop accuracy, using the results of refs. 11: 



m(m 2 ) 



= m I 



i+ C F a s (m 2 ) +{R 



2C< 



■>(t) ; 



(2.5) 



K(n f = 4) ~ 13.4; K(n f = 3) ~ 14.0. 
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3 Energy Shifts 



3.1 Order aj, A 4 /m 4 

Taking into account the expression for the Hamiltonian, Eq. (2.1), we can write 



E nl = 2m - + E 6 v E " 1 + 6 % E m + S npE nl . 



(3.1) 



Here the Sy^ E n \ may be easily evaluated with the formulas in the Appendix to rcf. 4. We define generally 
the analogue of the Bohr radius, 



mC'FO: s (p 
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and then, 
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(3 C F a 2 s (a 2 ) 



+ 



C F a s 21 + 1 -An 



m 2 n 4 (2Z + l)a 3: 
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log^+^(n + / + l) 



log^ + i>{n + l + 1) 



e(l) 



(3.2a) 
(3.2b) 
(3.2c) 

(3. 2d) 

m n 3 {2l + l)a 2 ' (3 ' 2c) 

We recall that constants are collected in the Appendix. For the vector states (T, T', Y"; J/tjj, tp', . . .) one 
has to add the hyperfine shift, at tree level, 



+^{n + l + l) 2 +i>'(n + l + l) 



Mdf , ^ 2r(n-l) rpi + 2 + j 

[ ] r{n + l + l) jKn -l-j-l) 
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(3.2f) 



The calculation of the second order contribution of V± L \ S^ L) E n i, is nontrivial, and may be found in the 



Appendix. We define 

S v , L) E nl = -m^^r [No +N, log — + 3 log — j- 
and then one has, for the lowest states, 



(3.3a) 
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(2,0) _ 1 - 2 7f 
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(3.3b) 
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432 12 



4 36 

- 3 - 



-CALCULATION OF QUARKONIUM SPECTRUM AND nib, m c TO ORDER a£- 

In addition to this one has to consider the nonperturbative (NP) energy splittings. The dominant ones 
are associated with the gluon condensate and are^ 



SnpE„i = me n in 2 Tr{a s G 2 ) (j-^j 



e n in & ir{a s G' 2 ) 
(mCi?a s ) 4 



(3.4) 
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9 945 ' 



Because (a s G 2 ) ~ A 4 , this is of order (A/m) 4 albeit with large coefficients: for all terms we have a fourth 
power of a s in the denominator, and for n > 1 the n 6 in the numerator of Eq. (3.4) grows very quickly out 
of hand. In fact it is the size of this term that limits the range of validity of our type of ab initio calculation. 

3.2 Higher corrections 

Besides the corrections reported in the previous subsection there are a few pieces of the higher order 
corrections that are known. First of all we have the relativistic, 0{v 2 ) corrections to the one loop potential W . 
These produce corrections of higher order, aP, but they are logically independent of the three loop ones that 
would produce terms of the same order but, presumably, smaller because of the extra 1/7T 2 characteristic 
of radiative corrections. These corrections may be incorporated and then can be considered to give an 
indication of the error committed in neglecting higher order perturbative corrections. They produce, for the 
ground state, the energy shift (ref. 4; typos corrected in ref. 8), 



lloop^^lO - M + A s , 
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(3.5) 




Next we have higher order NP corrections. These include finite size corrections, estimated in ref. 7, and 
contributions of higher dimensional operators, some of which were evaluated in ref. 9. The last produce the 
shifts, 

SnP, higher = ZIKTFi ~ \« M n ) 0^6 



m 5 (C F a s )» 

, .„ „, 141912051712 , ,„ , 484859657191424 
Ml,0) = — — — — , h(2,0) = 



844421875 



Oe(/i) 



h(2,l) = 



2040039729 
102150951135870976 

765014898375 ' 



(3.6) 



\qq\0) 2 



and may be used to estimate the size of the higher order NP contributions. For the quark condensate the 
vacuum saturation approximation is assumed, and the value of k is taken from ref. 10. For (G 3 ) one takes 
the value 0.065 GeV 6 . Anyway, these quantities are poorly known. 

It is important to realize that both (3.5), (3.6) should be taken as indications. With respect to the first, 
there is no guarantee that the coefficient of the three loop correction is not so large that it offsets the factors 
of 1/ir; indeed, this already happens to two loops where the coefficient is large, &i ~ 24. With respect to 
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(3.6), and apart from the fact that it does not include all the higher dimensional operators (those associated 
with size corrections are neglected 2 ), it is clear that one cannot consider rigorously a contribution 0(A e /m 6 ) 
so long as the radiative corrections to the 0(/l 4 /m 4 ) terms are not known. Nevertheless, we consider (3.5) 
and (3.6) as very useful for estimating the theoretical uncertainties of our calculation. 

4 Numerical Results 

Using the formulas deduced above one can evaluate the spectrum of heavy quarkonium systems. In principle 
one should take m, A, (a s G 2 ) from other sources and predict the masses of the quarkonium states. In practice 
it is better to use the known masses of the lowest states (T and J/tp) to evaluate the quark masses. The 
reason is that this produces by far the more precise evaluation available at present of these parameters, 
especially in the case of the b quark. The other parameters we take from independent sources. For the QCD 
parameter A we take, throughout this section, 

A(n f = 4, three loops) = 0.23±g;8§ GeV [ ct a (M%) ~ 0.114±g;85 ] . (4.1a) 

and for the gluon condensate, very poorly known, 

(a s G 2 ) = 0.06 ± 0.02 GeV 4 . (4.1b) 

Another matter to be discussed is the choice of the renormalization point, \x. As our equations (3.2, 3) 
show, a natural value for this parameter is 

,x=— , (4.1c) 

na 

for states with the principal quantum number n, and this will be our choice. For states with n = 1 the 
results of the calculation will turn out to depend very little on the value of /j,, provided it is reasonably close 
to (4.1c). Higher states are another matter; we will discuss our choices when we consider them. 

4.1 The 10 state of bb and the mass mi, 

As stated, we select, for the Y state, [i — 2ja. We then use Eqs. (3.1-4) to obtain the values of the b quark 
mass. To make apparent the contribution of the higher corrections, we have performed calculations taking 
into account only 0{a 2 s ), 0(a\), O(al) + \v\ x 1 loop terms, and finally the full 0(a 4 ) evaluation. The results 
are reported in Table I below 3 , where the errors correspond to the errors in Eqs. (4.1a, b). 





0{a 2 s ) 


O(ajj) lloop + rel.* 0(a 4 ) 


M 2 (GeV 2 ) 


3.233 


4.940 ~ 2.5 7.019 


m b (GeV) 


4.752 


4.858 4.939* 5.015t^ {A) T 0.005 «a s G 2 }) (vary /i 2 by 25%) 


m 6 (m 2 ) = 


4.209 


4.307 4.382* 4.453«;^ (A) =p 0.004 «a s G 2 }) ^ (vary /i 2 by 25%) 


* T1 


TABLE I: determinations of the 6 quark mass with increasing accuracy, 
lis value is taken from ref. 4, extrapolated to A = 0.23 GeV, (a s G 2 ) = 0.06 GeV 4 . 



In the estimate of the errors, the condition u. = 2/ a is maintained satisfied when varying A while for the 
error due to the variation of \i the other parameters are kept fixed (i.e., one no more has then fi = 2/a). The 
dependence of nib on /x should be taken as an indication of the theoretical uncertainty of our calculation. To 



2 The reason for doing so is that, at least nominally, the contribution of operators associated with the size, ((9G(0)) 2 ), 
is of higher order in a s than the ones considered in (3.6). See ref. 9 for details. 

3 We have arranged in Table I the results in terms of powers of a s ; we could have arranged them in increasing number 
of loops. Cf. Sect. 5 for this. 
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estimate other theoretical errors in our evaluation we proceed as follows. We either calculate m b including 
the full 0(v 2 ) corrections to one loop, Eqs. (3.5). Then we get, 

(1 - loop + rcl.) + (2 - loop, NR) : t = 7.009, m b = 5.010; m b {m 2 b ) = 4.448; (4.2a) 

or we may include the contribution of higher dimensional NP effects, as in (3.6). Then, 

with higher NP effects: m b = 5.018; m b {m 2 b ) = 4.455. (4.2b) 

We consider that the best result is that of O(af) reported in Table I, and take the difference with the 
quantities given in Eqs. (4.2) as a further estimation of the theoretical error of the calculation. In this way 
we get our best estimate, 

m b = 5.015±8;J^ (A) T 0.005 {{a s G 2 }) (vary fi 2 by 25%) ± 0.006 (other th. uncertainty) 

m b {m 2 b ) = 4.453±g;8?e ( A ) T 0.005 ((a s G 2 )) (vary ^ 2 by 25%) ± 0.005 (other th. uncertainty). ^'^ 

The values of a s (p, 2 ), a s (p 2 ) corresponding to fj, 2 = 7.019 GeV 2 are 

a s {n 2 ) = 0.24 ,a s (/i 2 ) = 0.40. 

The "theoretical" error coming from higher dimensional operators and higher order perturbative terms, 
Eqs. (4.2), is comfortably smaller that the errors due to the uncertainty on A, {a s G 2 ). We will henceforth 
omit these errors, so as not to double count them, and consider that the theoretical error is only that due 
to varying ^ 2 by 25%. If we now compose all the errors quadratically, then we obtain the estimate reported 
in the Introduction, Eqs. (1.1). 

4.2 M(T) - M{r\b); the decay Y -> e+e" 

The evaluation of refs. 4, 8 for the hyperfine splitting, and the decay of the T into e + e~ does not change, 
except that the favoured values of A, (a s G 2 ) and m b are now somewhat different. This improves slightly the 
agreement with experiment for the decay rate. The expressions are, 



M{V) Mir,) = m ^W^)^ 2 ) 3 (J + ^ + ^ 



x {1 + 



2 



(log^-l) +f (\ogC F & s + l) + B 



Os_ iigi iT{a s G 2 ) 
7T + 8704 m 4 56 



(4.4a) 



r(T - e+e") = r(°) x [1 + 5 wf + (5 NP ] 2 (1 + «5 rad ), 



r (0) = 2 



QbdQED 



2 

(mCF&sip 2 )) 5 (4.4b) 



M(Y) 

AC F a s 3/3 (, a[i \ a s 

Orad = ; <3wf = — j- I log — - 7E 5 

7T 4 V 2 / 7T 

r _ 1 [ 270 459 , 1838 781 1 K(a s G 2 ) 
° NP 2 L 108 800 " r 2 890 000 J m 4Q,6 - 

The corrections here are fairly large, particularly the radiative correction! 12 ! S ra d- Because of this the 
calculation is less reliable than what one would have expected. With the values of m b found in (4.3), 
one has the numerical results, 

M(T) - M( V ) = 48.4™ {A) ± 4 4 » ((a s G 2 )) t 7 w 2 2 (/i 2 = 7.019 ± 25%) (4.5) 

and 

r(T e+e-) = 1.135±g;5i {A) ± 0.11 ((a s G 2 )) ±g;^ {p? = 7.019 ± 25%). (4.6) 
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Note that, when varying A, (a s G 2 ), we have varied raj, according to Eq. (4.3), but we have not varied nib 
when varying /i. 

Higher order NP corrections due to the higher dimensional operators introduced in Eq. (3.6) are also 
known for the decay rate (see ref. 9). They read 

_ w(n) _ 1670626488940208128 

"NNP 2m 6 (CF<S) 10 485563688671875 ' 

Size corrections, however, are not known now. Snnp would produce a shift in the decay rate of ~ 0.11 keV, 
smaller than the contribution of (a s G 2 ) or the uncertainty caused by e.g. varying /x as in Eq. (4.3). We do 
not include Snnp either in the evaluation or the error estimate. 

The result for the decay is in reasonable agreement with experiment, 

r exp .( T -> e+e ~) = 1-320 ±0.04 keV. 
Composing the errors we obtain the figures quoted in the Introduction, Eqs. (1.2, 3). 

4.3 Higher states (n = 2) of bb 

The masses of the states with n — 2 will be next determined. As is clear from the expressions (3.2, 3) the 
natural choice of scale is now (j, = 1/a which gives n 2 = 3.05 GeV 2 . If we take this, adding or subtracting a 
25% to estimate the dependence of the calculation on the choice of scale then we obtain the results 

M (20, th) - M (20, exp) = 293^ McV (p 2 = 3 - 048 ± 25% )> 

M(21, th) - M(21, exp) = 174±^ MeV {p 2 = 3.048 ± 25%). ' ' ' 

We only present the errors that follow from variation of the scale /x 2 by 25%; slightly smaller ones are 
produced by the errors of A, (a s G 2 ). We do not explicit this: because of the size of the errors in Eq. (4.7) 
there is no point in going for a more detailed error analysis. 

Although they have decreased from th one loop evaluations (e. g. ref. 4), the errors are still fairly large 
here; within them, there is compatibility between theory and experiment. Agreement to a few MeV for both 
states is obtained if choosing /x 2 = 0.75/a ~ 2.3 GeV 2 or keeping /lx = 1/a and taking (a s G 2 ) = 0.036 GeV 4 : 
unlike for the states with n = 1 we have now strong dependence of the results on the parameters of the 
calculation. This is due to the large size of the corrections, perturbative and (especially) nonperturbative. 
This last is made more apparent when considering contributions of higher dimensional operators^ 9 !, which 
get completely out of hand for nl = 20 and are very large for nl = 21. In this context, it is satisfactory to 
realize that it is for this last state (21) for which agreement with experiment is best, and theoretical errors 
smaller. 

We will not discuss here the spin and tensor splittings among the states with nl = 21. The inclusion of 
three loop corrections to the potential only affects their calculation in that the preferred value for nib will 
be different now, which is a minute effect compared with the uncertainties of the calculation: one should 
realize^ 8 ] that, while the NP corrections to the energy levels with principal quantum number n contain 
a coefficient n 6 /a 4 , wave functions at the origin get a factor ~ n s /a®. This of course is what makes the 
calculation of M(T) — M(r)b) and the decay T — > e + e~ much less reliable than that of M(T) (or, equivalently, 
mf,) aud what makes the evaluation of tensor and spin splittings with n — 2 somewhat marginal. All one 
can do here is fit fi to the data; this is the procedure followed in ref. 8, and we have nothing new to report 
on this. 

4.4 The 10 state of cc and the mass m c 

The value of the parameter A used now, corresponding to that in Eq. (4.1a), is 

A(n f = 3, three loops) = 0.30±g;g§ GeV. 

The values for the c quark mass, deduced from the J/ip mass are then, 

0(a A s ) : t = 2.623, m c = 1.884, tti c (t77 2 ) = 1.547 
(1 - loop + rcl.) + (2 - loop, NR) : t = 2.611, m c = 1.875; tti c (t77 2 ) = 1.539. (4.8) 
with higher NP effects : t = 2.634, m c = 1.891; m c (m 2 c ) = 1.554. 
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Including errors we obtain the best estimate, analogous to that in (4.3) for the b quark: 




-0.089 
+0.145 



-0.096 
+0.156 



(varying u 2 by 25%) ± 0.011 (th. uncertainty) 
(varying /z 2 by 25%) ± 0.010 (th. uncertainty). 



(4.9) 



As is obvious from these equations, the errors are now much larger than for the b quark case, but our 
determination of m c still competes in accuracy with those based on QCD sum rules. 

5 Discussion 

The calculations of this paper are rather straightforward; but there are a few points that merit further 
discussion. First of all, our values for the quark masses are somewhat larger than existing estimates based 
on sum rules; for the MS masses, of 100 ~ 300 MeV, cf. ref. 10. In our opinion this is due to the influence 
of the terms of order a^, a 4 which we take into account, but which the sum rule evaluations, that stop at 
0(a 2 ), do not. Thus we consider our estimates to be the more precise and reliable ones. 

Secondly, and from Table I, it may appear that the series is diverging: from the 0(a 2 ) to the 0{a 3 s ) 
evaluation, mt increases by 106 MeV but from the last to the O(o 4 ) the increase is of 157. Actually, 
convergence is reasonably good. The increase between 0(a^) and 0(a 4 ) is due to two independent factors: 
inclusion of the two loop corrections to the potential, responsible for 46 MeV, and the relativistic corrections. 
Of these, 64 MeV for tree level corrections and 40 MeV for the mixed one loop-velocity correction. Each of 
the effects is small. Thus, if we included velocity corrections at every loop, the variation from zero to one 
to two loops would be diminishing. This is apparent if we compare the value obtained in ref. 4 (corrected 
for the increased values of A, (a s G 2 ) we are using now) with our results, with a variation of only 60 MeV. 
One can see this more clearly if we arrange the calculation in increasing number of loops including, at every 
step, the pertinent relativistic corrections 4 or including all loop corrections, but in increasing order of the 
velocity corrections: 



Finally, we remark the satisfactory stability we now have against variation of the renormalization scale, /z. 
This stability of the results against changes of /j, is made apparent by the fact that even multiplying or 
dividing the central value /i 2 = 7.019 by a factor of two only alters the central value of rrib = 5.015 GeV by 
98 MeV. The stability is due mostly to the inclusion of two loop effects; but attention should also be paid to 
the stabilizing influence of the NP corrections. These corrections are larger for larger n., exactly the opposite 
to what happens to the perturbative corrections. One could even fix optimal values of fx as those where the 
combined perturbative-NP effects would show a minimal dependence on /i. This is essentially the procedure 
adopted in ref. 4. Here, and because of our much smaller dependence on fi, we need not have recourse to 
such methods. 



4 Nevertheless, it is true that the two loop correction is large: fei ~ 24, bi/iv 2 ~ 2.5 and a 2 s b\/-K 2 ~ 0.14: to be 
compared with ai/ir ~ 0.5, a s a\ji: ~ 0.11. We are clearly near the limit of convergence of the perturbative series. 
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Appendix 
Constants 

O = U- |n /; /3i = 102- fn/ 

a _ 2847 5033 „ , 325^,2 

31C A - 20T F n f C F - 2C A 

di = ~ 1.47; a 2 = — -2.33; 

36 2 

_ [me + » f(3) ] - [f - 16 C (3)] C F T F n f + *§T F n}} 

~ 24.30; 

14C F - 21C A -16C F + AC A 
as = § , a 4 = , 

a 5 = 2C F + ^-^+47>log2 ; 
3 15 

B = |(1 - log2)I> - lT F n f + 11Ca ~ 9 ° F ~ 0.29. 

Second order contribution 

We now calculate d^.( L) E n i. A simple, and rather accurate approximation may be obtained with use of the 
following trick. The second order shift given by a potential W may be written as 

6$E nl = (nl\W\S^R nl ), 

where 



\6WR nl ) = P nl \R%); 



P n i is the projector orthogonal to the state nl and 



\R { nh = M> + E tp \ \n'l)(n'l\W\nl) 

77 - &n' 



is the wave function to first order. The trick is to use for this not the result of a Rayleigh-Schrodinger 
formula, but that obtained from a variational principle. For our case (cf. ref. 4, Eqs. (79-81)), 



R W- 2 / ("-*- 1)' (*\ /nb 2l+1 , , , 



log(na/x/2) + ^(n + Z + l)-l 
1 H (hots 



6- 1 = a- 1 

For the ground state, this gives, after a trivial calculation, 

c(2) j-, PlC F ala s r 2 a M 1 2 

V^ 10 = ' m 16^ i l0g T _ 7E • 

This simple method gives correctly the coefficients of log a^i, log 2 a[i and misses the constant term by ~ 10%. 
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-CALCULATION OF QUARKONIUM SPECTRUM AND nib, m c TO ORDER CKg- 

For the exact calculation one uses the representation of the Coulombic Green's function given e.g. by 
Voloshin 5 . We write 

\(n,i\vl L) \k,i)\ 2 



c(2) 



r (L) Enl 



E 



E n — Ek 



and the sum over k includes an integral over the continuous part of the spectrum. Instead of doing this 
computation directly we have used the more general function 



V \{n,l\vl L) \k,l)f _ \(n,l\v} L) \n,l)\ 2 v \(n,l\v} L) \k,l)\' 

k k^n 



E - E k 



E - E n 



+ A {0) +A (1) (E-E n ) + ... 



(2) 

and then 5 v{L) E n i = ^4(o)- For the function defined above one uses the representation of the Coulombic 

Green's function given e.g. by Voloshin in the second article in rcf. 6 [Note that there is a misprint in 
formula (15) there, and (s + I + 1)! must be changed to (s + 21 + 1)!]. In this way we get for the different 
N n i defined in Eq. (3.3), and for arbitrary quantum numbers, 



(n-l-l)\ 



N 

1-1-2 

E 



(n,l) 



iP(l + l + n)-l 



(s + 2/ + l)! 



(n + l)\ 



E 



■n)[V(l 



0-2] 



(n + l)\ ^ s!(s + Z + l-n) 3 (n - I - 1)\ ^ (s + 2 1 + 1)! (s + I + 1 - n) 3 ' 

s=0 v / \ / s — n —l v ' K ' _ 



5 An estimate neglecting the continuum contribution was been given in ref. 14. The approximation is not good, 
although since the quantity is small its effect in the evaluation of m b was not important. 
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